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The Nature of the 10 Kilosecond X-ray Flare in Sgr A* 
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Abstract. The X-ray mission Chandra has observed a dramatic X-ray flare - a brightening by a factor of 50 for 
only three hours - from Sgr A*, the Galactic Center supermassive black hole. Sgr A* has never shown variability 
of this amplitude in the radio and we therefore argue that a jump of this order in the accretion rate does not 
seem the likely cause. Based on our model for jet-dominated emission in the quiescent state of Sgr A*, we suggest 
that the flare is a consequence of extra electron heating near the black hole. This can either lead to direct heating 
of thermal electrons to T c ~ 6 • 10 11 K and significantly increased synchrotron-self Compton emission, or result 
from non-thermal particle acceleration with increased synchrotron radiation and electron Lorentz factors up to 
7c il 10 5 . While the former scenario is currently favored by the data, simultaneous VLBI, submm, mid-infrared 
and X-ray observations should ultimately be able to distinguish between the two cases. 
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accretion disks - black hole physics 



1. Introduction 

Sgr A*, the compact radio core at the center of our 
Galaxy (Reid et al. 1999; Backer & Sramek 1999), has 



been perplexing modelers since its discovery (Balick & 
Brown 197^ ). In contrast to nearby LLAGN jHo 1999| ) 



Sgr A* was until recently only positively detected as a ra- 

~|6; 



- 0.01 pc ( 


Hallcr et al. 1996; 


Eckart & Genzel 1996 




Ghcz 


|et al. 1998 


) and its integrated radio luminosity has re- 



mained steady within a factor of two (Zhao, Bower, & 
Goss 2001), at ~ 10~ 9 orders of magnitude less than its 
corresponding Eddington luminosity. All models to ex- 
plain the radio emission so far have focused on radiative in- 
efficiency as the primary explanation for this dimness, and 



is best fit with a power-law (spectral index a ~ 0.3), 
which is significantly harder than that of the quiescent 
state (a ~ 1.2). The longest time scale (10 ks) corresponds 
to ~ 390r s where r s = 2GAf,/c 2 is the Schwarzschild ra- 
dius, which argues against thermal bremsstrahlung from 
the outer radii, e.g. from a standard Advection Dominated 
Accretion Flow (ADAF; [Narayan et al. 1998| ). The small- 
est time scale in the flare is roughly 600 s, suggesting ac- 
tivity at scales of ~ 20r s , which means the flare originated 
close to the central engine. 

The variability and the spectral index of Sgr A* in 
the X-rays are consistent with synchrotron self-Compton 
(SSC) from the innermost regions near the black hole, e.g., 
the nozzle of a jet (FM00; |Yuan, Markoff, fc Falcke 2001[) 



are comprised mainly of accretion/inflow solutions ( Melia 



or a magnetic dynamo within the circularized accreting 



Liu, fc Coker 2001 ; Narayan et al. 1998 ) outflow solutions 
( Falcke, Mannheim , & Biermann 1993c Falcke & Markoff 



plasma (Melia, Liu, & Coker 2001). In this picture, the X 



200C, hereafter FM00) a nd combinations thereof ( Yuan 



Markoff, fc Falcke 200l|). A recen t review of Sgr A* can 
be found in [Melia fc Falcke (200l| ). 



Recently, Sgr A* was finally det ected in the X-rays 
by Chandra ( Baganoff et al. 2001b| ) with a r ather soft 
spectrum. Du ring the second observational cycle, Baganoff 
et al. (2001a ) detected an X-ray flare lasting about 10 
ks and with a peak luminosity ~ 50 times higher than 
the quiescent state ( Baganoff et al. 20011: ). The averaged 
flare spectrum after taking into account dust scattering 



rays are inverse Compton up-scatter ed synchrotron pho- 
tons from the so-ca lled submm-bump ( Serabyn et al. 1997 ; 
Falcke et al. 1998 ). Since the submm-bump is thought 
to be produced close to the black hole, very short time 
scale variability (several hundred seconds) was already 
predicted (FM00). In the following we would like to ex- 
plore the various scenarios which could lead to a dramatic 
X-ray flare within the jet model. 



2. Models 
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We start with our basic jet emission model (|Falcke fc 
Biermann 1999| ; FM00), consisting of a conical jet with 
pressure gradient and nozzle. The parameters in the noz- 
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zle for the quiescent state are determined from the un- 



derlying accretion disk as described in Yuan, Markoff, & 



Falcke (2001), and as summarized below. All quantities 
further out in the jet are solved for using conservation of 
mass and energy, and the Euler equation for the acceler- 
ating velocity field. We take the distance to the Galactic 
center as d gc = 8.0 kpc. 



observations. This argument does not hold for the poorly 
sampled data at other wavelengths and we only consider 
single-epoch measurements which most likely only reflect 
the quiescent Sgr A* spectrum. 

The effects of the M-flare and the T e -flare can be 
modeled simply by changing the jet power and tempera- 



CleHi'ly, in order Lo produce an X-ray flare, one or sev- 
eral parameters had Lo have suddenly changed in Sgr A + . 

In Figs. 2 and 3 in FM00, we showed how the radio and 
X-ray spectra in the jet model change if one changes the 
magnetic field - a similar result would be expected for 
an increase in particle density - or the electron temper- 
ature by a small amount. The former would be expected 
for an increased jet power or accretion rate, which would 
result in simultaneous flaring at all frequencies with lit- 
tle change in spectral index. In the latter scenario, how- 
ever, the X-rays flare much stronger with a hardening of 
the spectrum, because SSC is very sensitive to changes in 
electron energies. This type of fast heating could in prin- 
ciple occur via instantaneous transfer of energy from the 
magnetized plasma in the accretion flow to the radiating 
particles, e.g. as would be expected from the sudden dis- 
charge of energy in magnetic flares through reconnection 
(e.g., [Biskamp 1997| ). 

On the other hand, we know that non-thermal parti- 
cle distributions are quite common in jets in AGN and 
X-ray Binaries (XRBs), leading to the appearance of op- 
tically thin power laws in the spectra. Observations of 
jets in XRBs (e.g., |Fender 200 lty and some AGN (e.g., 



ture, respectively, in our published models (FM00, Yuan 



Meisenheimer, Yates, & Roeser 1997) seem to hint at a 
common type of power law with typical spectral index of 
a ~ 0.6 — 0.8. While the exact mechanism is not yet firmly 
established, and reconnection may also contribute, first or- 
der diffusive shock acceleration leads more naturally to an 
electron distribution with the index p ~ 2 — 2.6 depe nding 
on the shock co mpression ratio (4j| oc E~ 



see e.g., lones 



& Ellis an 1991 ) . Such accelerated particles would result in 
a significant increase of optically thin synchrotron emis- 
sion, with spectral slope a = (p — l)/2. 

In the following we therefore explore three 



scenarios 



for thd origin of the X-ray flare: increased jet power or 

accretion rate, increased heating of relativistic particles, 
or sudden shock acceleration of the particles. We will refer 
to these three models as the M-flare, the T e -flare and the 
shock-flare, respectively. 

3. Results 



Markoff, & Falcke 2001). We assume that the jet carries 
away a fixed fraction of the accretion energy Mc 2 , and 
that this energy is divided evenly between the kinetic en- 
ergy carried by the cold plasma, and the internal energy 
carried by the magnetic field and hot electrons. Once the 
electron temperature T e in the nozzle is fixed, assuming a 
Maxwellian distribution, the jet nozzle density no is de- 
termined via approximate equipartition from the magnetic 
field Bq. In the quiescent state, the relevant parameters 
for our most recent fit are T c =~ 2 • 10 11 K, n ~ 9 • 10 5 
cm -3 and B ^20 G. Fig. |l| shows the prediction for a) 
the M-flare, with jet power (oc B 2 ) raised by ~ 3 via in- 
creasing the jet nozzle magnetic field Bq to ~ 35 G while 
holding T e fixed (which in turn increases no by ~ 3 to 
~ 3 • 10 6 cm -3 ) and b) the T c -flare for T e raised by a fac- 
tor of ~ 3 to T c ~ 6 • 10 11 K, while holding no and Bq 
fixed. The parameters were chosen to match the amplitude 
of the X-ray flare data, shown with its error box as well. 
For comparison we show in the figure also the quiescent 
jet+disk spectrum. 

As expected, the M -flare strongly over-predicts the ra- 
dio flux by a large factor. In fact, such a huge flare in the 
radio has never been reported and in addition, the spec- 
tral index is far too steep. The T e -flare fares much better: 
the predicted radio flux is close to already observed ra- 
dio flare maxima and the X-ray spectrum becomes very 
hard during the flare. The model also predicts significant 
brightening in the MIR range during the X-ray flare event, 
due to the shift of the submm-bump to higher frequencies, 
which should exceed currently available non-simultaneous 
MIR/NIR limits. In contrast to the radio, the MIR regime 
has not been sampled well enough to decide whether such 
flares exist. However, Genzel fc Eckart (1999 ) and Scrabyn 
et al. (1997) report observations where Sgr A* could have 



been detected during a brief period with unusually high 
flux densities at 350 /im and 2.2 fim. Clearly, this needs 
to be confirmed and reassessed in light of the new X-ray 
observations. 

The shock-flare scenario requires more discussion, as 
it involves the effects of diffusive shock acceleration in the 
jet. This has been done already by Markoff, Falcke, & 



Since no simultaneous radio or mid-infrared 



(MIR)observations are available we include in our 



Fender (2001), where the scaled version of the jet model 
previously used to explain Sgr A* (Yuan, Markoff, & 



figures an "upper radio envelope", showing the highest Falcke 2001; FM00) has successfully been applied to X 



flux ever detected at each radio freque ncy in long-term 
monitoring of Sgr A* with the VLA (Zhao, Bower, fc 



Goss 200 1| ). While it is possible that this type of X-ray 



flare is so rare that it was never before captured by radio 
observations, it seems statistically unlikely given the huge 
radio database compared to only two cycles of Chandra 



ray binaries in the low/hard state by including shock ac- 
celeration. Because the low/hard state is characterized by 
a very faint, possibly AD AF- like accretion disk as in Sgr 
A*, the ambient photon field is not strong enough to re- 
sult in significant inverse Compton (IC) cooling, allowing 
shock accelerated electrons to achieve rather high energies. 
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Fig. 1. Fit of the jet model (solid line) to the flare data of Baganoff et al. (2001a) (a) with the M-flare, raising the 
jet power by a factor of ~ 3 , and (b) with the T c -flare by raising the temperature of the electrons by a factor of 3, 
compared to the quiescent jet and disk model (dashed line). The radio data and IR upper limits are from the data 
set compiled and presented in Melia & Falckc 2001, where IR data are from single-epoch observations only. The upper 
radio points show the highest flux detected at that particular frequency, as compiled by Zhao, Bower, & Goss (2001). 
The lower X-ray data show the quiescent state spectrum of Baganoff et al. (2001b). 



Following [Markoff, Falcke, fc Fender (200l[ ) the par- 
ticles would be accelerated up to a maximum energy 
-Ee,max = 7e,mai"»eC 2 , which is reached when the syn- 
chrotron loss rate equals that of acceleration. We use the 
simple parallel shock acceleration rate 



acc j 



eB 



Ushy 



(1) 



where u s h is the shock speed in the plasma frame. The 
parameter £ < c/3 e /u sn ( Jokipii 1987 ) is the ratio between 
the parallel diffusive scattering mean free path and the 
gyroradius of the particle, and ranges from a lower limit 
at £ = 1 up to typically a few 10 2 (e.g., Jokipii 198^ ). For 
a magnetic field of ~ 20 G as found in our model of the 
quiescent state, the acceleration time scale is ~ 0.1 sec 
for even j c = 10 5 electrons, and hence is shorter than the 
dynamical time scale at the black hole. 

Setting the standard synchrotron loss rate = 

|cT7c/3e^r = ^acci we can solve for the maximum elec- 
tron energy achieved by acceleration 7 c . max . If we define 
as a reference value £ = £2 100, the maximum synchrotron 
frequency is then 



0.29i/ c ~ 1.2 • 10 



20t-l f u sh\' 
42 \~) 



Hz 



(2) 



where v c ~ j^7 2 max (eB) / (m c c) is the critical synchrotron 
frequency. This value is not dependent on the magnetic 
field, the jet power, or the shock location as long as we 
are in the synchrotron cooling dominated regime. 

Because the shock accelerated particles responsible for 
the X-ray synchrotron will have very high energies (70^ 
10 5 ) for the low magnetic fields further out in the jet, the 
synchrotron cooling time scale will be very short, on the 
order of ~ 10 2 s. This means that re-acceleration along the 



jet is required to maintain the population, and will result 
in rapid cooling if the acceleration is switched off. We thus 
approximate the shock acceleration as continuous starting 
at a distance z sn . For X-ray binaries we found that the 
shock acceleration must begin relatively close to the nozzle 
at z sh ~ 10 ~ 10 2 r s ([Markoff, Falcke, fc Fender 2001 



Markoff et al., in prep.). This location is determined from 
the data by extrapolating the synchrotron X-ray curve to 
where it meets the optically-thick, flatfish spectrum in the 
radio-IR. This intersection is unique for a fixed spectral 
index, and gives z sn because the self-absorption frequency 
scales inversely with z in the jet model. 

If we then fix for simplicity the fraction of accelerated 
part icles at 50% and keep the other p arameters as in FM00 
and |Yuan, Markoff, fc Falcke (200l[ ) we can calculate the 
resultant shock-flare model spectrum as shown in Fig. ^| 
As the spectral index becomes harder for a fixed X-ray 
flux, the optically thick turnover must occur at lower fre- 
quencies, i.e. further out in the jet. For a standard spectral 
index of a ~ 0.8 as typically seen in AGN, the shock accel- 
eration region must be at ~ 16r s , which is consistent with 
the observed time scales. However, the assumed standard 
AGN spectral index is only marginally compatible with 
the spectrum observed for the X-ray flare, which poses a 
problem for such a model. Taking on the other hand the 
reported best-fit X-ray spectral index at face value would 
imply a = 0.3 and require z sn ~ 10 4 r s . This is very far 
in comparison to other jet systems and furthermore ruled 
out by the observed short time scales. 



4. Discussion 

We are able to explain the 10 ks flare in Sgr A* detected by 
Chandra by heating the radiating electrons within the jet 
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Fig. 2. Fit to the flare data of [Baganoff et al. (2001a| ) for 
the shock- flare model, other data the same as Fig. 1. 



model of FMOO, either so they remain quasi-thermalizcd 
(T c -flare) or via the non-thermal process of shock acceler- 
ation (shock- flare). A flare due to an increase in accretion 
rate (M-flare) appears very unlikely because of the gener- 
ally low level of radio variability. 

Of the two remaining scenarios, the shock-flare is more 
intriguing because it offers a solution where the X-ray flare 
occurs without a great effect on the lower-frequency emis- 
sion, consistent with the observed lower radio variability 
of Sgr A* over the last decades. At the marginal end of 
the fit, it can also explain the shortest variations via the 
location of the shock or fast radiative cooling, and predicts 
a spectral index consistent with that seen in other AGN 
systems. Although the radio flux does not change much 
for this case, the presence of the optically thin tail would 
predict a significantly larger radio profile (more extended, 
optically thin jet emission; see FMOO for a discussion of 
this point) on the sky and a shift of the centroid of the 
radio emission. However, in the radio astrometric work of 
Reid et al. (1999 ) no such shift has been detected so far. 

Alternatively, the T e -flare with its sudden heating of 
hot (T c ~ 6 • 10 11 K) electrons by, e.g., magnetic recon- 
nection, can explain the X-ray flare via increased SSC 
emission, similar to models for the quiescent state spec- 
trum. The fast variability can be explained by the small 
source size and outflow with v ~ c, leading to fast adi- 
abatic cooling, while radiative cooling is not as impor- 
tant (t syn ~ 5 • 10 4 for T c = 6 • 10 11 K electrons in the 
submm-bump). In contrast to the shock-flare model, the 
Te-flare model fits the reported X-ray spectrum much bet- 
ter. However, the radio variability is larger than in the 
synchrotron case, but still falls along the "envelope" of 
highest radio fluxes observed so far (Fig. |l|). In addition 
the model predicts simultaneous MIR flaring, in a regime 
where no monitoring data is currently available. So al- 
though the T -flare case is favored over the shock-flare 
case in terms of the fit to the X-ray flare data, only simul- 



taneous submm/MIR/X-ray and VLBI observations in the 
near future will unequivocally determine its viability. 

For the T e -flare, assuming the protons have at least 
the same temperature one can compare the energy den- 
sity of the plasma and the gravitational binding energy, 
GM,m p n/ R, yielding 



TnkhT 



GM,m p nR 



3.4 



GoOsec) 



T 



6 • 10 n K 



(3) 



for a M, = 2.6T0 6 M Q black hole and a relativistic plasma 
with r = 4/3 (e.g., Konigl 1980) at a distance R = ct 
from the black hole set by the variability time scale. Under 
such simple assumptions the plasma would not be grav- 
itationally bound, which is certainly consistent with the 
jet scenario. 
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